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Abstract: To plan and design safe and reliable space missions, it is necessary to take into account 
the effects of the space radiation environment. The environment during large solar energetic 
particle events poses the greatest challenge to missions. As a starting point for planning and 
design, a reference environment must be specified representing the most challenging 
environment to be encountered during the mission at some confidence level. The engineering 
challenge is then to find plans and mission design solutions that insure safe and reliable 
operations in this reference environment. This paper describes progress toward developing a 
model that provides such reference space radiation environments at user-specified confidence 
levels. 
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Abstract: To plan and design safe and reliable space missions, it is 
necessary to take into account the effects of the space radiation 
environment. The environment during large solar energetic particle 
events poses the greatest challenge to missions. As a starting point 
for planning and design, a reference environment must be specified 
representing the most challenging environment to be encountered 
during the mission at some confidence level. The engineering 
challenge is then to develop plans and mission design solutions that 
insure safe and reliable operations in this reference environment. 
This poster describes progress toward developing a model that 
provides mission-specific reference space radiation environments at 
user-specified confidence levels. The examples shown here are for 
episode-integrated proton fluence spectra. 


Step 1: Select solar energetic particle events for the data 
base 

• NOAA Criterion 

o Onset: The first of 3 consecutive data points with >10 MeV 

proton fluxes > 10 PFU (protons/cm2.ster.sec). 

o End: The last data point >10 PFU*. 
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Identifying events and finding onset and end times. 


October 20, 1994 
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The smallest events are affected most by the NOAA criterion. 



Data 

Sources: 



Step 2: Determining the energy spectra for each event 

• Determine the spectra in the instrument-defined energy bins 

• Convert the spectra to a common set of energy channels. 

o requires spectral fitting to redistribute the energy into the 
common set of energy bins 
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We have fit 214 episodes with: 

- The Ellison-Ramaty Model (90) 

- The Band Function (116) 

- The Weibull Function (136) 

In many cases two or three of these models gave good 
fits. 






Step 3: Construct cumulative distributions of the 
fluxes for each energy channel. 

• Below is an example of the cumulative distribution for the 
21.89 MeV energy channel 

• Each cumulative distribution must be fit by a Frechet 
Distribution 

p (*> = (C - f b )/(C - £) 



Step 4: Convolve the cumulative distributions with a 
Poisson distribution to create an extreme value 
distribution for the fluence in each energy channel. 

F t (S = exp{-f<r[l - p(0]} 

• The event frequency varies over the solar cycle so extreme value 
distributions must be constructed for each year of the solar cycle. 



Summary 

* We have developed a model for estimating the worst-case 
episode-integrated proton spectrum that is: 

o Specific to the mission start date and duration 

o At a user-specified confidence level 

* Limitations : Can’t obtain high confidence level reference 
environments for long missions that are supported by data. 

* We plan to extend this model to: 

o Alpha particles and Heavy ions 


Average Number of SPEs per Year 

(based on 5 solar cycles) 


to 


18 

16 

14 

12 




-6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 


Year relative to Maximum 


Step 5: Construct the Worst Case Spectrum 

• Using the user-supplied launch date and mission duration, determine 
from the figure above, the average number of SPEs expected during 
the mission. 

• Use the mission-average number of SPEs to construct the extreme 
value distribution for each energy bin. 

• Use the user-supplied confidence level to determine the fluence in 
each energy bin that corresponds to this confidence level. 

• Construct the mission-specific worst-case fluence spectrum. 

• An example is shown in the next column. 

o It is for a three-year mission launched in 20 1 6 at a 90% confidence 
level. 


o We plan to construct similar models for peak fluxes and 
mission-integrated fluences. 
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Determining a reference worst-case 
environment. 

Selecting the SPEs. 

Determining the energy spectra for 
each event. 

Constructing the cumulative spectra. 

Finding the extreme-value distributions. 
Finding the worst-case energy spectra. 



What is a reference environment 
and why do we need one? 


• A Reference Worst-Case Environment in the local 
interplanetary medium is needed for: 

- Spacecraft design 
— Mission planning 

• Objective: Determine an environment that: 

— Won't be exceeded at a specified confidence level 

- For a user-specified mission, e.g. 

• Launch date, mission duration, heliospheric location 

• This talk focuses on the proton component of the 
Reference Worst-Case Environment 



Past work on reference 
environments 

• JPL Model Objective: 

• Xapsos Models 

• Rosenqvist Model 

• Jiggens Model 




How is the Reference Worst-Case 
Environment determined? 



• Build a database of SPEs with measured 
elemental spectra and time-histories. 

• Decide how to treat solar cycle dependence of 
SPE. 

• Form cumulative frequency distributions for each 
differential energy channel. 

• Construct extreme value distributions for each 
channel. 

• Using the user-supplied orbit, launch date 
mission duration and confidence level, construct 
the worst case elemental spectra. 



The first step: Event Selection 



• Use a consistent selection criterion that can be applied over 
the broadest time span => Use GOES data (1974-present) 

• Group sequences of events that appear to be correlated into 
episodes (same active region, sympathetic triggering, etc.) 

• Minimize biases due to how onset and ending times are 
determined. 

We have chosen to use the NOAA criterion: 

• Onset: The first of 3 consecutive data points with >10 MeV 
proton fluxes > 10 PFU*. 

• End: The last data point > 10 PFU*. 


*protons/cm 2 .ster.sec. 


Counts/cm 2 sec sr 


Applying the NOAA Criterion, example 1 



• Missing or bad data: 6:00-6:15 on 2/6 

• Multiple events: There appears to be an episode of >2 correlated events 

• 1 or 2 Events: Strictly applying the criterion, the first event ends at 14:00 
on 2/6, but this is not really the end of the episode. 



Counts/cm 2 sec sr 



8 8.1 8.2 8.3 8.4 8.5 8.6 8.7 8.8 

Day of November 

• Where is the onset? Strictly, at 7:50 but The event started much earlier 
at ~3:25. 


Counts/cm 2 sec sr 



Day of February 

• Where does it end? Strictly, at 16:20 but the trend falls below 10 PFU at 
about 17:30. 


Counts/cm 2 .ster.sec 
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The smallest events are affected 
most by the NOAA criterion. 




Episode-integrated proton spectra 


Proton Data Sources 

- 7/74-10/01: IMP-8 GME 

• 167 episodes 

- 11/02 - Present: GOES 

• 47 episodes 

Common spectral binning 

- Adopted 29 GME channels 

Redistribute GOES data into 
GME channels 

- Fit GOES spectra with a spectral 
form, e.g. double power law 

- Distribute flux into GME 
channels 





Data sources for other 
elemental spectra. 


* Helium 

- 7/1974 - 10/2001: IMP-8 GME 

- 11/2001 - present: GOES 

• Heavy Ions 

- 7/1974 - 10/2001: IMP-8 CRNE 

- 11/2001 - present: ACE SIS 



Putting all the proton spectra in 

a standard format 


• The IMP-8 GME format with 29 energy bins 
was chosen as the standard format. 

— The spectra for the events before Nov. 2001 are 
taken from GME measurements. 

• For events after October 2001, GOES data 
were used. 

- The 7-energy-bin GOES spectra were fitted. 

- The best spectral fits were used to re-bin the data 
into the GME format with 29 bins 



Measured proton spectra and fits 
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• We have fit 214 episodes with: 

— The Ellison-Ramaty Model (90) 

— The Band Function (116) 

— The Weibull Function (136) 

• In many cases two or three of 
these models gave good fits. 



Cumulative spectra 


• Cumulative spectra were constructed for each energy 
bin. 

— An example is shown below. 





Solar Cycle variation is SPE 

probability 


• The probability of occurrence of solar particle events 
varies over the 11-year solar cycle as shown below. 


Average Number of SPEs per Year 
(based on 5 solar cycles) 
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Extreme Value Model 


• The number average number of SPEs, N, expected during 
the mission is found by summing over the part of the solar 
cycle covered by the mission. 

• N is used to construct a Poisson distribution, P(n). 

Pin ) = [e { ~ N \N) n ]/n\ 


The cumulative distributions, C((/>), are fit with the 
following function (Xapsos et al., 1998): 

C(0 = if b - f b )/{f b ~ f b ) 
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These distributions are convolved to get extreme value 
distributions for each energy bin: 

Fi(</>) = exp^NP t (</))} 


Protons/cm 2 .ster.MeV 


Worst-Case Proton Spectrum 



Energy in MeV 


Summary 

We have developed a model for estimating the 
worst-case episode-integrated proton spectrum 
that is: 

— Specific to the mission start date and duration 

— At a user-specified confidence level 

Limitations : Can't obtain high confidence level 
reference environments for long missions that are 
supported by data. 

We plan to extend this model to: 

— Alpha particles and Heavy ions 

We plan to construct similar models for peak 
fluxes and mission-integrated fluences. 


